Analysis and identification of ferroptosis-related diagnostic markers in rheumatoid arthritis

Abstract Background Rheumatoid arthritis (RA) is an autoimmune, inflammatory joint disease. There is growing evidence that ferroptosis is involved in the pathogenesis of RA. This study aimed to search for diagnostic markers of ferroptosis in RA and to analyse the potential mechanisms and clinical value. Materials and methods RA-associated datasets were used from the publicly available GEO database. Three methods of machine learning were applied to screen biomarkers. The diagnostic efficacy of the results was also verified by receiver operating characteristic (ROC) curve, external dataset, qRT-PCR and Western blot. Enrichment analysis was performed in this process, while protein–protein interaction (PPI) analysis and immune infiltration correlation analysis were performed using biomarkers, and competing endogenous RNA (ceRNA) networks were constructed to search for prospective therapeutic targets. Results MMP13 and GABARAPL1 can be used as ferroptosis diagnostic genes in RA. The ROC curve and PPI result demonstrated that MMP13 and GABARAPL1 had an excellent diagnostic value. The results of signature genes in the external dataset, qRT-PCR and Western blot further confirm our findings. The enrichment analysis showed that p53, MAPK and NOD-like receptor signalling pathways may be involved in the process of ferroptosis in RA. In addition, two ferroptosis diagnostic genes in RA participate in the occurrence of ferroptosis in RA via oxidative stress, metabolism and immune response. Immune infiltration analysis showed that RA extensively infiltrated B cells, T cells, macrophages and other immune cells. Persistent immune activation may be an essential reason for the progression of ferroptosis in RA. We also obtained five potential therapeutic agents for RA and some long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) regulating ferroptosis diagnostic genes. Conclusions Our study suggests that MMP13 and GABARAPL1, which are closely linked with oxidative stress and immunological modulation, can be used as ferroptosis-related potential diagnostic markers in RA and provide new clues regarding the diagnostic and therapeutic targets of ferroptosis in RA.


Introduction
Rheumatoid arthritis (Ra) is an autoimmune disease that attacks multiple joints throughout the body and can cause irreversible joint damage and disability [1].swelling and inflammation of the joint synovium, destruction of cartilage, production of autoimmune antibodies, and systemic multi-organ lesions are the main pathological features of Ra [2].Nonsteroidal anti-inflammatory medications, disease-modifying antirheumatic drugs and glucocorticoids are among the medications used, which can reduce the symptoms or slow the course of Ra. however, their complex pathogenic mechanisms and adverse drug reactions vastly limit the implementation of treatment options [3].these medications cannot alter the nature of Ra; they can only partially ease the patient's symptoms.current evidence suggests that ferroptosis in fibroblast-like synoviocytes (Flss) and synovial giant cells is associated with Ra development, and ferroptosis also triggers chondrocyte destruction in Ra [4,5].therefore, finding the diagnostic markers that precisely regulate ferroptosis in Ra is of great clinical importance in guiding the therapeutic strategy for Ra.
Mitochondrial structure and function abnormalities are the primary morphological hallmarks of ferroptosis, a regulated cell death [6]. the pathomechanism is manifested by increased reactive oxygen species (ROs) due to lipid peroxidation and massive intracellular iron deposition, ultimately leading to cellular dysfunction [7].Many studies have shown that releasing intra-articular proinflammatory factors in Ra patients promotes iron deposition in the synovium, maintaining a persistent synovial inflammatory state [8].ROs and the severity of Ra are correlated favourably. in recent years, the development of high-throughput sequencing has allowed us to recognize the development process of disease from multiple dimensions, such as somatic mutations and altered genomic expression [9,10].it was found that specific ferroptosis-related genes (FRGs) are involved in the progression of Ra.For example, Xiang et al. found that the FRG, slc2a3, had elevated expression in Ra tissues by analysis and identification, and showed down-regulation in response to ferroptosis inducers (Rsl3), thus causing the occurrence of ferroptosis in Ra [11].as an FRG, RRM2 has also been found to be differentially expressed between Ra patients and normal subjects, and the viability of damaged cells is also increased with the inhibition of this gene, it is thought that FoxO signalling pathway and inherited metabolic disorders play a role in this [12].VeGFa, PtGs2 and JUN are also thought to be involved in the regulation of ferroptosis in Ra by interacting with immune cells [13].this implies that a few critical genes play a crucial role in the ferroptosis of Ra.Ferroptosis has received increasing attention in the field of Ra research in recent years.exploring the relationship between ferroptosis and Ra at the level of gene targets and identifying corresponding therapeutic strategies seem to be an up-and-coming area.some studies have explored ferroptosis diagnostic markers in Ra, but fewer studies have validated their diagnostic value.Based on bioinformatics, ferroptosis signature genes (FsGs) of Ra were identified by machine learning and experimental validation methods. in this study, we first obtained the differentially expressed FsGs between Ra and normal samples by three datasets, Gse29746, Gse55235 and Gse55457, and validated their diagnostic value using the external dataset Gse77298.then, we further validated the results obtained based on the GeO dataset by qRt-PcR and Western blot in vitro experiments.the study will help explore the potential biological targets associated with ferroptosis in Ra, analyse their regulatory mechanisms, and provide a different strategy for treating Ra.

Data acquisition and pre-processing
Four Ra-related mRNa transcriptome datasets (Gse29746, Gse55235, Gse55457 and Gse77298) were downloaded from the GeO (https://www.ncbi.nlm.nih.gov/geo/) database.the dataset was processed using the R package as follows: Gse29746, Gse55235 and Gse5545 were corrected using the 'limma' and 'sVa' R packages to remove batch effects [14].then, log2 processing is performed on the data with large values, and finally, the data are merged.all the combined data were analysed with relative logarithmic expression (Rle) using R software, and the results were visualized with boxplot (supplementary Figure 1).the training group included 32 Ra synovial samples and 30 normal synovial samples. in addition, Gse77298 were used as independent sequences as the validation group, including Ra/normal synovium samples, were 16/7, respectively.eight hundred and forty-four FRGs were downloaded from the ferrdb database (http:// www.zhounan.org/ferrdb).

Acquisition of ferroptosis differentially expressed genes (DEGs)
the expression of FRGs was extracted from the training group using the 'limma' R package.DeGs for ferroptosis between Ra and normal synovial samples were then extracted using the 'limma' R package as defined by p values <.05, and the results were visualized using volcano plot and heat maps.correlation analysis of DeGs was performed using the 'corrplot' R package.

Enrichment analysis
correlation analysis of DeGs' biological functions and signalling pathways was conducted using the cluster-Profiler 'R' package [15].the process used Gene Ontology (GO) and the Kyoto encyclopedia of Genes and Genomes (KeGG).GO consists of three components, biological processes (BPs), cellular components (ccs) and molecular functions (MFs).p < .05 was used as the setting condition.

Candidate FSGs selection
We used three machine learning algorithms to screen for reliable FsGs in Ra. lassO regression analysis was performed to select variables by a regularized analysis algorithm to identify genes that significantly differed between Ra synovium and normal synovium.the sVM-RFe algorithm selects and visualizes the most relevant FsGs through a non-linear kernel and accurately interprets the direction and strength of the association between the FsGs and the results [16].WGcNa can find the most relevant gene modules for the disease by a weighted approach [17].the 'glmnet, ' 'e1071R, ' and 'WGcNa' R packages were separately used for analysis to recognize the genes that are of diagnostic significance.the genes obtained from the three methods were crossed to get FsGs for subsequent analysis and validation.

Confirming the reliability of RA FSGs for diagnosis and gene set enrichment analysis (GSEA)
after obtaining the FsGs crossed by the three algorithms, the receiver operating characteristic (ROc) curves of the FsGs were plotted using the pROc R package, and their diagnostic value was assessed based on the area size of the ROc curves.Meanwhile, to verify the accuracy of these FsGs as ferroptosis diagnostic genes for Ra, we used an independent dataset, Gse77298, to validate the expression differences of the FsGs.Finally, we performed Gsea analysis on the FsGs.

Construction of PPI network of FSGs
to explore the intrinsic association of FsGs with FRGs and pathways, a PPi network with medium confidence (minimum required interaction score of 0.4) was obtained using the stRiNG database to elucidate the potential mechanism of action.

Immune infiltration correlation analysis
the differences in the relative abundance of 22 infiltrated immune cells in normal and Ra synovial samples from the training group were analysed by the ciBeRsORt algorithm [18].spearman's correlation analysis was performed on the infiltration of immune cells and FsGs.Vioplot and ggplot2 R packages plot violin plots and correlation heat maps to analyse and visualize the results.

Prediction of potential drugs for RA and construction of competing endogenous RNA (ceRNA) network
Based on FsGs, we identified potential therapeutic agents for Ra by using the gene-drug interaction database (DGidb) (https://www.dgidb.org/)[19].to explore the interactions between FsGs and post-transcriptional regulators long non-coding RNa (lncRNa) and microRNa (miRNa), we constructed a ceRNa network between the three using miRanda, miRDB, targetscan and spongescan bioinformatics tools.the regulatory network was analysed and visualized using cytoscape_v3.9.1 software [20].

Extraction and cultivation of synovial cells
in this study, five pairs of Ra and normal synovial tissue were collected from the First affiliated hospital of Guangxi Medical University (Nanning, china) between June 2023 and september 2023.synovial tissue for Ra was obtained from five patients who fully consented to knee arthroplasty and met the diagnostic criteria for acR [21].the control group was from five patients who required arthroscopic surgical treatment for anterior cruciate ligament injury or meniscal injury.the synovial membrane was taken from the normal area within the joint.all patients were hospitalized in the First affiliated hospital of Guangxi Medical University.Written informed consent of each patient was obtained during the sample extraction process.Detailed clinical details are provided in supplementary table 1.Briefly, synovial tissue was rinsed 2-3 times with phosphate-buffered saline (PBs) and then well minced and digested in DMeM medium (Gibco, Pleasantville, NY) containing type i collagenase (1 mg/ml; Gibco, Pleasantville, NY) for 4-6 h.Next, synovial cells were filtered using a strainer and centrifuged to obtain synovial cells.the cells were cultured in a DMeM medium containing FBs (20%, Gibco, Pleasantville, NY) and penicillin/ streptomycin (1%, solarbio, Beijing, china) in an incubator at 37 °c with 5% cO 2 .the experimental cells were third-generation synovial cells.all procedures of this experiment were approved by the ethics committee of the First affiliated hospital of Guangxi Medical University (approval number: 2023-e206-01).

RNA extraction and qRT-PCR
RNa was extracted from Ra and normal samples according to the instructions of the RNeasy™ Plus animal RNa extraction Kit (Beyotime, R0032, shanghai, china), and 1 μl of the sample was tested for the concentration and purity of the extracted RNa by UV spectrophotometry.a total of 1 mg of RNa was then reverse transcribed to cDNa in a total volume of 20 μl by a reverse transcription tool (takara, Beijing, china).Real-time qPcR was then performed according to the specifications of sYBR-Green (thermo Fisher scientific, Waltham, Ma) operating instructions.the 2 −ΔΔct method was used to analyse, and the gene GaPDh was used as a standardized measure.independent replicate experiments were performed three times, and each set of samples was repeated thrice.Primer sequences are shown in table 1.

Western blot
Protein was extracted from synovial cells using a configured protein lysis buffer, including phosphatase inhibitor cocktail (comWin Biotech co., ltd., Beijing, china), protease inhibitor (Medchemexpress, Monmouth Junction, NJ), RiPa lysate (solarbio, Beijing, china) and PMsF (solarbio, Beijing, china).after the loading buffer was obtained, the protein concentration was determined by Bca protein assay (Beyotime, shanghai, china).Protein samples were separated by 7.5% and 15% polyacrylamide gels and transferred to polyvinylidene fluoride (PVDF) membranes, then closed with protein-free fast blocking buffer (shanghai, china) for 30 min.the membranes were soaked in the corresponding primary antibody GaPDh (Proteintech, shanghai, china), MMP13 (Proteintech, shanghai, china), GaBaRaPl1 (abcam, cambridge, UK) and incubated overnight in a shaker at 4 °c.Western blot was visualized using Proteinsimple's Wes system.Gray scale values of Western blot were analysed using image-J software (BioRad, hercules, ca).

Statistical analysis
all bioinformatics analyses were performed on R v. 4.1.13.Differences between the two groups were compared using the Wilcoxon test.correlations between different variables were analysed using spearman's correlation analysis.experimental data from in vitro experiments were analysed using the unpaired two-tailed t-test method.a statistically significant difference was considered at p < .05.

DEGs associated with ferroptosis and enrichment analysis
We obtained the combined cohort (Gse29746, Gse55235 and Gse5545) after eliminating batch effects using the 'limma' and 'sVa' R packages.From 844 FRGs downloaded through the ferrdb database, we identified 389 FRGs in the sequence set.We further screened 144 DeGs related to ferroptosis in Ra in the training group, of which 60 were upregulated and 84 were down-regulated.Figure 1(a,B) (volcano plot and heat map) shows the distribution of DeGs among different samples.We used KeGG and GO enrichment analysis to demonstrate the biological functions and signalling pathways associated with the DeGs.KeGG showed the relationship between DeGs and MaPK/NOD-like receptor (NlR) signalling pathway, p53 signalling pathway, ferroptosis and apoptosis (Figure 2(a)).elevated expression of P53 inhibits slc7a1 1/GPX4 expression and thus induces cellular ferroptosis [22].Ferroptosis is a cell death regulated by intracellular iron deposition and lipid peroxidation and is closely related to the intracellular antioxidant system. in GO (Figure 2(B)), DeGs were significantly enriched in cellular response to oxidative stress, ROs metabolic process, positive regulation of cellular catabolic process, and other BPs, indicating that DeGs are closely associated with inflammatory response processes related to iron and lipid metabolism. in cc, DeGs were enriched in autophagosome, NaDPh oxidase complex, cytoplasmic vesicle lumen and secretory granule lumen. in MF, DeGs were associated with oxidoreductase activity, superoxidegenerating NaD(P)h oxidase activity, iron/ferrous binding and protein serine/threonine kinase activity, suggesting that DeGs are involved in processes of cellular iron metabolism and antioxidant regulation.the key to ferroptosis lies in intracellular iron deposition and the weakening of lipid antioxidant function.enhancing the activity of antioxidant system and inhibiting iron deposition are effective measures to treat ferroptosis in Ra [23].

Screening and identification of FSGs
We first screened 389 FRGs in the training group using the WGcNa analysis method.Different modules were identified according to the correlation between  similarly, we screened out 20 DeGs using the sVM-RFe algorithm (Figure 3(D)).Finally, the genes obtained from the above three methods were crossed to get 12 FsGs (slc1a5, GaBaRaPl1, eGFR, MaPK8, PRKaa2, il1B, sNca, aDaM23, PtGs2, BeX1, tRiB2 and MMP13) (Figure 3(e)).Following that, we constructed ROc curves to assess the diagnostic value of these FsGs and screened them according to the criteria of excellence defined by aUc value ≥0.85 [24].the results showed that GaBaRaPl1, eGFR, tRiB2 and MMP13 met the inclusion criteria, with aUc values of 0.946, 0.901, 0.858 and 0.857, respectively (Figure 4(a)).We next validated the 4 FsGs using the external dataset Gse77298 (p < .05 was deemed significant).the results showed that only two FsGs, MMP13 and GaBaRaPl1, passed the validation and were consistent with the expression trend of the training group (Figure 4(c-F)).We also constructed model ROc curves for these two FsGs with an aUc value of 0.959 (Figure 4(B)).We also examined the diagnostic accuracy of these two FsGs in the Gse77298 cohorts.the aUc values of MMP13 and GaBaRaPl1 were 0.946 and 0.768, respectively (Figure 4(G,h)).the above results showed that these two FsGs provided good diagnostic accuracy.

Co-expression network of FSGs and ferroptosisrelated proteins
to elucidate the potential relationship between the FsGs and ferroptosis-related proteins, we constructed a PPi network diagram between two FsGs, ferroptosis-related proteins (Fth1, GPX4, slc7a11 and P53) and related pathways through the stRiNG database (Figure 5).We found a close interaction between the FsGs and ferroptosis-related proteins.

Gene set enrichment analysis
Based on the expression profile of the FsGs, Gsea analysis revealed the functional and pathway enrichment of the FsGs in the high and low expression groups.We found that the up-regulation of MMP13 expression was closely related to the cytokine receptor interaction, graft versus host disease, leishmania infection, adaptive immune response, B cell activation, immune effector process, leukocyte mediated immunity and regulation of cell activation (Figure 6(a,B)).this suggests that upregulated MMP13 may mediate ferroptosis in Ra by activating the body's immune response.More interestingly, the down-regulation of GaBaRaPl1 was significantly enriched in chemokine signalling pathway, adaptive immune response, meiotic cell cycle, nuclear chromosome segregation, side of membrane, chemokine activity and cXcR chemokine receptor binding (Figure 6(c,D)).this suggests that FsGs play a role in immune response, cell cycle and metabolism, cytokine interactions, protein synthesis and metabolism.all of these are associated with ferroptosis in Ra. analyses of the expression of FsGs have improved our understanding of the pathogenesis related to ferroptosis in Ra.

Immune cell infiltration
an autoimmune disease affecting joints throughout the body; Ra is distinguished by the infiltration of immune cells and the production of autoimmune antibodies [25].We found differences in B-cell memory, plasma cells, macrophages M1 and mast cells activated by analysing 22 immune cell infiltrations between Ra and normal samples (Figure 7(a)).the expression of most immune cells is upregulated in Ra.We also plotted the heat map between FsGs and immune cells.the results showed that FsGs were positively correlated with macrophages M0, t cells cD4 memory activated, dendritic cells resting and NK cells activated, and negatively correlated with dendritic cells resting and B cells naive (Figure 7(B)).

Identification of the expression levels of two FSGs
in this study, two FsGs were detected by qRt-PcR and Western blot to identify their expression between Ra patients and normal subjects.compared with normal samples, the mRNa expression of MMP13 increased in Ra samples, while GaBaRaPl1 decreased (Figure 9(a)).Western blot measures protein levels of MMP13 and GaBaRaPl1 in synovial cells.When comparing Ra samples to normal samples, higher protein levels of MMP13 and lower protein levels of GaBaRaPl1 were found in the former (Figure 9(B)).

Discussion
Ra is a long-lasting, inflammatory joint condition.the invasion of synovial inflammation and cartilage destruction is often accompanied by adverse joint results and ultimately leads to disability.as a disease with high activity, the production of autoantibodies and the appearance of early joint damage imply the importance of early diagnosis and intervention.therefore, the search for useful diagnostic markers and new therapeutic approaches is clinically important for early recognition and intervention of the disease.Many studies have found that ferroptosis, a programmed cell death distinct from apoptosis and characterized by iron deposition and lipid peroxidation, plays a vital role in the pathological process of Ra [26,27].however, the relationship between FRGs and Ra and the mechanism of action has not been clarified.in the study, two FsGs of Ra were identified.the findings of immune infiltration analysis and the interactions between lncRNa, miRNa and FsGs have triggered our thinking on the relationship between Ra and ferroptosis.Ra is currently a disease with no established etiology and pathogenesis.however, pertinent research has demonstrated that alterations in specific pathways and persistent immunological activation are intimately linked to the incidence and progression of Ra [28].enrichment analysis showed that 144 ferroptosis DeGs were significantly enriched in MaPK/NlR/p53 signalling pathway and involved in iron deposition and oxidation regulation processes.the enrichment result suggests these genes play an important role in Ra ferroptosis.activation of the MaPK inflammatory pathway has been shown to down-regulate the expression of slc7a11 and GPX4, leading to ferroptosis [29].Xie et al. found that activation of MaPK pathway contributes to ferroptosis of cancer cells [30].Meanwhile, li et al. also found that the ferroptosis pathway mediated by ROs/MaPK and p53 can kill nasopharyngeal carcinoma cells [31].as a crucial regulator of ferroptosis, the dysfunction of p53 in Flss is associated with synovitis and the destruction of bone and cartilage.sun et al. indicated that elevated expression of p53 could induce ferroptosis in cells through the p53/slc7a1/GPX4 pathway [22].however, there are few studies on the correlation between the expression of p53 and ferroptosis in Ra [32].as a family of cytoplasmic proteins, the NlRs family plays a crucial role in pathogen recognition and innate immune response.By identifying the appropriate ligand, the NlR signalling pathway can activate the downstream signalling pathway NF-κB and then promote the transcription of the inflammasome NlRP3.Next, NlRP3 induces the activation of caspase-1, which promotes the processing and secretion of precursors of il-1β and il-18 [33].a crucial ferroptosis regulator is GPX4. in Ra, the overactivation of NlRP3 may limit the production of GPX4, resulting in a significant buildup of ROs and iron-dependent cell death [34].several pathways, including MaPK and apoptosis, are stimulated by the NlR signalling pathway at the same time [35].the activation of a series of pathways and the release of inflammatory factors initiate the immune response process.We built a PPi network to investigate the connections between FsGs, ferroptosis-related proteins and associated pathways.the findings suggest that FsGs may have a role in ferroptosis directly or indirectly via related pathways, but the precise mechanism remains to be further investigated and discussed.
MMP13 is the protease that plays a major role in Ra and is produced by chondrocytes [36].synovitis of Ra releases inflammatory cytokines that stimulate MMPs production.MMPs can degrade collagen and aggrecan of the extracellular matrix of the cartilage, thereby causing irreversible destruction of the articular cartilage, bone and tendon.cartilage fragmentation in turn irritates synovium, leading to aggravation of inflammation [37].the buildup of ROc also can activate MMP13.the elevated expression of MMP13 can cause more severe joint damage in Ra patients [38].By interacting with the NF-B binding site on the MMP gene, the transcription factor aP-1 can increase the expression of MMP [39].this suggests that MMP13 may be related to the NlR signalling pathway and involved in the process of ferroptosis.however, the mechanism of MMP13 in ferroptosis in Ra is unclear and was screened as FsG in this study, which piqued our curiosity.GaBaRaPl1 belongs to one of the subfamilies of the yeast autophagy-related eight families, which has the function of transporting vesicles and proteins, regulating energy metabolism and oxygen radical homeostasis, and is involved in various pathological processes such as autophagy, cell proliferation and death [40]. it has promise as a molecular target for autoimmune or inflammatory diseases [41].autophagy, particularly selective ferritinophagy, increases lipid peroxidation and iron accumulation, resulting in cellular ferroptosis [42].Down-regulation of GaBaRaPl1 in hepatocellular carcinoma cells diminished their sensitivity to erastin, a ferroptosis inducer, according to Du et al. [43].interestingly, down-regulated GaBaRaPl1 expression was also found in whole blood of Ra patients [44].We believe that GaBaRaPl1-mediated autophagy can cause cells to undergo ferroptosis.the role of these two FsGs for ferroptosis in Ra remained unclear and was identified for the first time as diagnostic markers of ferroptosis in our study.the diagnostic value of these two FsGs was also demonstrated by ROc analysis.this was also confirmed in our in vitro experiments, in which MMP13 showed significant upregulation of expression in Ra samples, while GaBaRaPl1 showed significant downregulation of expression.the above results suggest that these two FsGs can be used as potential ferroptosis-related diagnostic markers for Ra.
Ra is an autoimmune disease in which the infiltrated immune cells in the joint and the existing autoantibodies participate in the whole process of the disease [45]. in the study, Ra samples were shown to have higher B cells memory, plasma cells and macrophages M1 infiltration.FsGs are associated with macrophages M0, t cells cD4 memory activated, dendritic cells resting, NK cells activated and B cells naive.according to studies, the leading causes of joint damage are MMP-1 and MMP-13 released by synovial cells in Ra [46].at the same time, it has been shown in Ra that lymphocytes and monocytes promote the expression of MMPs and cytokines utilizing microparticles [47].according to the available data, the GaBaRaPl1 protein is crucial for adaptive immunity because it processes and presents antigens via Mhc-i and Mhc-ii molecules, promoting the production of cD8+ or cD4+ t cells [48].additionally, the Mhc-i protein in dendritic cells uses the autophagy of GaBaRaPl1 to cross-present antigens [49].activation of immune cells such as t cells, B cells, dendritic cells and macrophages promotes the production of rheumatoid factor and acPa, leading to increased synovial inflammation and irreversible joint destruction [50].the uncontrolled adaptive immune response and activation of innate immunity will exacerbate the progression of ferroptosis in Ra.Various immune cells, mainly B cells, t cells and macrophages, play an essential role in the ferroptosis development of Ra, and our study yielded consistent results [45].We also identified five potential therapeutic agents, including apratastat, prinomastat, doxycycline hyclate, doxycycline and doxycycline calcium.Moreover, to determine which lncRNas compete with FsGs for miRNa binding, we obtained an interrelationship network of two mRNas, 53 lncRNas and 47 miRNas.Previous studies have indicated that has-miR-19a can reduce the severity of Ra through negative regulation of MMP13 [51]. the negative regulation of GaBaRaPl1 by has-miR-19a and has-miR-145 affected the proliferation and differentiation of mesenchymal stem cells and endothelial progenitor cells [52,53].these findings indicate the potential for new ferroptosis therapeutic targets in Ra. therefore, we believe that MMP13 and GaBaRaPl1 engage in the process of ferroptosis in Ra and can serve as valuable diagnostic markers.Besides, they have a close relationship with pathways related to ferroptosis with further potential as therapeutic targets.however, more studies are necessary to explain the interaction between diagnostic markers and ferroptosis.admittedly, we have to acknowledge the shortcomings of this study.the sample size for analysis and verification is small in the first place.second, there is a need for further studies to explore the detailed mechanism of the diagnostic markers in Ra ferroptosis.

Conclusions
through bioinformatics and in vitro experiments, two FsGs, MMP13 and GaBaRaPl1, were confirmed to be differentially expressed between Ra and normal patients and could serve as valuable diagnostic biomarkers.Persistent immune activation and excessive oxidative stress as well as the regulation of ferroptosis-related pathways may be the possible mechanisms of ferroptosis caused by diagnostic biomarkers in Ra. the lncRNas and miRNas associated with diagnostic markers also provide new molecular targets for the treatment of Ra.
different gene-traits (Figure3(a)).the absolute value of correlation was taken and the module gene with the highest absolute value had the highest correlation with Ra. therefore, the Meturquoise gene module, which includes 215 ferroptosis-associated genes, was considered (Figure3(B)).Using the lassO logistic regression algorithm, we screened 21 DeGs from the 144 DeGs related to ferroptosis in Ra (Figure3(c)).

Figure 1 .
Figure 1.ferroptosis deGs between RA and normal samples.(A) Volcano plot of upregulated and down-regulated genes.(B) Heat maps of deGs, red for up-regulation, blue for down-regulation and white for no difference.

Figure 3 .
Figure 3. screening of RA fsGs.(A) dendrogram of WGcnA gene modules.The top and bottom corresponding parts are the hierarchical clustering map of genes and the gene module map.(B) Module-trait correlation heat map.each row represents a module and each column represents a trait.The middle two columns represent traits of RA synovium and normal synovium.Blue is a negative correlation, orange is a positive correlation, and the darker the colour, the stronger the correlation.The p values of correlation and significance are labelled in the cells.(c) fsGs screened from deGs using the lAsso regression algorithm.(d) fsGs screened from deGs using the sVM-Rfe algorithm.(e) Venn diagram of intersecting genes obtained by the three methods (WGcnA, lAsso and sVM-Rfe).

Figure 4 .
Figure 4. data set validation and diagnostic value analysis of RA fsGs.(A) Roc curves to validate the diagnostic accuracy of RA's fsGs.(B) Roc curves of the fsGs (GABARAPl1, MMP13) model.(c, d) expression of the fsGs in the merged dataset.(e, f) expression of fsGs in the Gse77298 dataset.(G, H) Roc curves of the fsGs in the Gse77298 validation group.

Figure 5 .
Figure 5. exploring the relationship between fsGs and potential signalling pathways of RA or fRGs.

Figure 7 .
Figure 7. Analysis of immune cell infiltration.(A) The difference of 22 kinds of immune cell infiltration between RA group and normal sample group.(B) correlation heat maps between 2 fsGs and 22 types of immune cells.Blue means negative correlation, red means positive correlation (*p < .05;**p < .01;***p < .001).

Figure 8 .
Figure 8. Gene networks of RA potential drugs and ceRnA regulatory networks.(A) The potential drug-fsGs network of RA was constructed by cytoscape, with red representing genes and blue representing drugs.(B) ceRnA network regulation map of lnRnA-miRnA-fsGs.

Table 1 .
Primer information on relevant genes in the study.